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A Simultaneous Magneto-Dielectric Phase Transition in RbCoBrs 

Takashi Shirahata * and Tota Nakamura 

Department of Applied Physics, Tohoku University, Aoha, Sendai, Miyagi 980-8579 

(Received February 2, 2008) 

We have modeled magneto-dielectric phase transitions in ABXa-type layered triangular 
lattice compomids. The model consists of a spin system (magnetic transition) and a lattice 
system (dielectric transition) . Magnetic exchange interactions are supposed to change with a 
relative position between two spins. This assumption produces an effective coupling between 
the lattice and the spin. Applying the nonequilibrium relaxation method, we have found 
that a simultaneous magneto-dielectric phase transition occurs when energy scales of the 
spin part and the lattice part coincide. An intermediate phase like the partial disordered 
phase disappears. Both systems relax frustration of the other system cooperatively and 
realize ordered states. 

KEYWORDS: Magneto-dielectric transition, layered triangular lattice, nonequilibrium relax- 
ation method, relaxation of frustration by lattice distortion 



1. Introduction 

The ABXs-type layered triangular lattice antiferromagnets have been attracting great 
interests both experimentahy and theoretically. ^"^^^ As the typical compounds we may list 
CsNiCls, CsCoCls, RbMnBra, KNiCls and RbCoBrs. The crystal structure is hexagonal close 
packed. Face-sharing BXg octahedra run along the c-axis forming a BX3 chain. Magnetic B^+ 
ions form an equilateral triangular lattice on the c-plane. It causes frustration in the exchange 
interactions. Exchange interactions along the BX3 chains are much stronger than those on the 
plane. This spin system can be considered as a quasi-one-dimensional system with frustration 
on the c-plane. 

Successive magnetic phase transitions occur in the Ising antiferromagnet on the two- 
dimensional triangular lattice with ferromagnetic next-nearest-neighbor interactions.^' ^'^^ A 
low-temperature magnetic structure is the ferrimagnetic state (we abbreviate it with "Ferri" 
hereafter). There exists a partially-disordered (PD) phase between the paramagnetic (Para) 
phase and the ferrimagnetic phase. One of three sublattices is completely disordered in this 
phase. The other sublattices take antiferromagnetic configurations. The PD phase is also 
considered to exist in the layered three-dimensional model.^'^'^'^'^''^ Since no anomaly of the 
specific heat was observed in the experiments^ between the ferrimagnetic phase and the PD 
phase, there is a possibility that this change of the magnetic structure is a crossover between 
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Fig. 1. Typical crystal structures in ABXa-type compounds. Solid circles depict magnetic ions 
and open circles depict X~ ions, (a) A symmetric structure at high temperatures. The space 
group is P6c/mmc. We call this structure "lattice Para" in this paper, (b) A room-temperature 
KNiCls structure. The space group is PQscm. We call this structure "lattice Ferri" . (c) Another 
low temperature structure. The space group is P3cl. We call this structure "lattice PD". 

the PD-like state and the ferrimagnetic-like state. 

A typical lattice structure of the ABX3-type compounds at high temperatures is shown in 
Fig. 1 (a). The space group is P6s/mmc. The structure is sometimes referred to as a CsNiCla- 
type. Magnetic ions forming an equilateral triangular lattice sit on a level plane (c-plane). We 
call this structure "lattice Para" in this paper because it is a symmetric structure at high 
temperatures like the paramagnetic state. When the temperature is lowered, structural phase 
transitions occur. Each BX3 chain shifts upward or downward keeping the relative distance 
between atoms in a chain unchanged. Some compounds take a structure in which two of three 
sublattices on the triangular lattice shift upward while the other one shifts downward as shown 
in Fig. 1 (b). The space group is P63cm. It is a room temperature structure of KNiCla. Since 
the displacement is t-T-i type, we refer to this structure as "lattice Ferri" in this paper. Other 
compounds take a structure with one sublattice shifting upward, one shifting downward and 
the rest unchanged as shown in Fig. 1 (c). The space group is P3cl. We refer to this structure 
as "lattice PD" since the displacement is t-i-0- 

It is known that each BX3 chain possesses negative charge. When a structural phase 
transition occurs, it can be observed by dielectric polarization. For example, a structure of 
the room-temperature KNiCl3 family (Fig. 1 (b)) induces macroscopic polarization. This is a 
distinct difference compared with the other structures in Fig. 1. Morishita et al.^^^ classified 
structural phase transitions in ABX3-type compounds by the polarization observation. A 
structural phase transition and a magnetic phase transition occur at different temperatures. 
A structural phase transition usually occurs at a higher temperature. This is because an 
energy scale of the lattice system is larger than that of the spin system. A lattice changes its 
structure from the symmetric one to a distorted one. A magnetic state remains paramagnetic. 
When the temperature is lowered, a magnetic phase transition occurs from the paramagnetic 
state to the PD state or to the ferrimagnetic state. 
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Recently, Morishita et al. found that a structural phase transition and a magnetic phase 
transition occur at the same temperature in RbCoBrs.^"^^ The magnetic susceptibility shows 
a peak at 37K, which suggests that a magnetically ordered state is realized at lower tem- 
peratures. The magnetic state has not been identified. The structural phase transition is 
observed at the same temperature by a peak of the dielectric constant. The low temperature 
phase exhibits finite spontaneous polarization observed by a pyroelectricity measurement and 
a D-E hysteresis measurement. A possible lattice structure with finite polarization is the 
room-temperature KNiCls structure shown in Fig. 1 (b). Therefore, it is conjectured that the 
lattice changes its structure from the CsNiCla-type to the room-temperature KNiCls-type at 
this temperature. 

In this paper, we make it clear why both magnetic and structural phase transitions occur 
simultaneously. We propose a model Hamiltonian which explains this phenomenon in §2. The 
model consists of a magnetic part and a lattice part with an effective coupling between them. 
We perform nonequilibrium relaxation (NER) analyses on this Hamiltonian. A brief review 
on this method is given in §3. Numerical results are presented in §4. Summary and discussions 
are given in §5. We obtain a phase diagram with regard to a temperature versus a ratio of 
energy scales between a magnetic part and a lattice part. A possible scenario which explains 
the simultaneous magneto-dielectric transition is given. 



We construct a model Hamiltonian in order to explain the magneto-dielectric phase tran- 
sition in RbCoBr3. First, we define state variables. One lattice point possesses a spin variable 
Sij and a lattice variable aij. Here, the subscript i denotes a site in the c-axis, while j denotes 
a site on the c-plane. Since the magnetic ion Co has a strong uniaxial anisotropy, we treat it 
with an Ising spin as Sij = ±1. For the lattice system, we define 0"^ by a displacement from 
the symmetric lattice point realized in the high-temperature phase. It is known experimen- 
tally that Co ions move only along the c-axis. Therefore, we suppose for simplicity that aij 
takes +1, —1 or depending on whether an ion shifts upward, downward or remains at the 
symmetric point. There are six different states at each lattice point. 

We introduce the following Hamiltonian consisting of a lattice Hamiltonian and a spin 
Hamiltonian. 



A lattice Hamiltonian is supposed to take an expression of elastic energy with regard to the 
lattice variables. A spring constant is denoted by , 



2. Model 



— Tilattice + "^spin 



(1) 
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n.n. 

i (jk) 

Here, J^, j}' and are spring constants of the nearest- neighbor pairs in the c-axis, the 
nearest-neighbor pairs on the c-plane and the next-nearest-neighbor pairs on the c-plane, 
respectively, {jk) denotes a pair in the nearest neighbor and (jl) denotes a pair in the next- 
nearest neighbor. It is known experimentally^^) that lattice dimerization within a CoBra-chain 
is negligible compared with the whole chain shift. The lattice is considered to be very hard 
only along the c-axis. This evidence can be modeled by setting take a positive value and 
larger than the other spring constants. The lattice system possesses a quasi-one-dimensional 
aspect. We have also set < and J^^ > in order to realize an t-T-i conformation (lattice 
Ferri) at low temperatures. This choice is consistent with an exclusion volume effect. The 
following parameters are used in our simulations: 

J,L = 5,jL = -l,jL = 0.1. (3) 

The temperature is scaled by | |. 

A spin Hamiltonian is defined by the following expression with the exchange integrals 
dependent on the lattice variables. 

n.n. 

- y^y^l-^l - ^jfWij — (^ik\)SijSik 
i (jk) 
n.n.n. 

-EE (-^2 - A J| IcTi,- - au \)SijSa (4) 

, Jf and Jf are exchange integrals in the high-temperature symmetric lattice structure. 
We have supposed that deformation of the lattice always decreases magnitudes of exchange 
integrals. Since spins interact via X ions by the super-exchange mechanism, overlap of elec- 
tron orbitals determines strength of the interaction. The overlap decreases when an relative 
angle along the exchange path has changed even though the direct distance becomes short. 
Therefore, it is natural to consider that magnitudes of the exchange integrals decrease when 
the lattice is deformed. For simplicity, the decrease is supposed to be proportional to an ab- 
solute value of the difference of lattice variables: \uij — (Ti'j'\, and a coefficient is AJ^^ ^ 2)- The 
nearest-neighbor interactions within the c-plane are supposed to be antiferromagnetic and the 
next-nearest-neighbor interactions are supposed to be ferromagnetic in order to realize the 
ferrimagnetic state in the ground state. The interactions along the c-axis are set ferromagnetic 
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and stronger than those within the c-plane. We define a ratio of jf to jj" by k and use the 
following parameters in our simulations: 

= 5k, jf = -K, jI = 0.1k, {k = Jf/J\) (5) 

A4,i,2) = 0.2JS,i,2)- (6) 
Experimental estimates for the exchange integrals are Jf ~ 62K, Jf ~ 2.5K, and 
jf — IK.^^^ A value of is estimated by a position of a broad maximum peak of x\\- 
jf and jf are estimated using a molecular field approximation on the layered equilateral 
triangular lattice antiferromagnet,'^^ where the successive transitions of Para^ PD Ferri 
occur. An application of the theory is not trivial because a direct transition of Para — > Ferri 
occurs in the present system. The obtained results may have differences from true values. 
Therefore, we consider that our choices of the exchange integrals are not inconsistent with 
the experimental situation except for jf being relatively small. Monte Carlo simulations gen- 
erally become difficult when a value of jf is large. A lattice length along the c-axis should 
be enlarged in accordance with a value of jf . It requires much longer CPU time to perform 
Monte Carlo simulations. We have chosen a moderate value of jf = 5k in order to avoid such 
time-consuming situations. It is also noted that a value of A J^^ ^ 2) is very large: exchange 
integrals may be 20%-40% decreased by a lattice distortion. This parameter controls an effec- 
tive coupling between the spin and the lattice. When it is large, a lattice distortion strongly 
influences a spin structure through changes in the exchange integrals. The spin and the lattice 
is independent when AJ^^ =0. We have set the parameter large so that we can easily 
observe the spin-lattice cooperative phenomenon. We consider that an essential part of the 
present system can be drawn out by these choices of parameters. 

A parameter k is controllable and we have carried out simulations at 

K = Jf/J\ = 0, 0.5, 1.0, 1.3, 1.6, 2.0, 2.5, 00. (7) 

When K is small, the lattice system dominates the whole system. The structural phase transi- 
tion is considered to occur at a high temperature. If the room-temperature KNiCl3 structure 
(lattice Ferri) is realized, the spin order favors the PD state because one interaction bond is 
stronger than the other two bonds in a unit triangle as shown in Fig. 2(b). On the other hand, 
if the lattice structure of the P3cl symmetry (lattice PD) is realized, the ferrimagnetic order 
is favored as shown in Fig. 2(c). In either case, the ground state is the ferrimagnetic state 
because of finite ferromagnetic next-nearest-neighbor interactions. 

3. Method 

We observe the following two physical quantities. One is the sublattice order parameter: 

iVsub 



ml = -^Y^Y^ai, (8) 
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Fig. 2. (a) A fully-frustrated triangular lattice, (b) A distorted triangular lattice realized in the room- 
temperature KNiCls-family of Fig. 1 (b). Since the interaction between A and B is weaker than 
that between B and B, two B-spins form antifcrromagnetic state, while an A-spin takes an up or 
down state randomly. This is the PD configuration, (c) Another distorted lattice realized in Fig. 1 
(c). Since \Jab\ > \Jbb\, A- and B- spins form antiferromagnetic state. This is the ferrimagnetic 
configuration. 



i jerj 

where 77 = a, /3, 7 denotes one of three sublattices in the triangular lattice, and Ng^ih denotes 
a number of sites in one sublattice. m]^ is sublattice polarization and is sublattice mag- 
netization. The other physical quantity is the 1/3 structure factor defined using sublattice 
polarization/magnetization: 

/v3 = ^^EK)'-K"^i'+"^X")/2^ (10) 
/1/3 = ^(^E(^')'-(^>?'+"^>?")/2^ (11) 

These structure factors take definite values when the Ferri or the PD state is realized and 
vanish when the symmetric/paramagnetic state is realized at high temperatures, f^^ detects 
a phase boundary between an symmetric lattice phase and an deformed lattice phase, ff^^ 
detects a phase boundary between a magnetic order phase and the paramagnetic phase. 

The sublattice order parameters distinguish the Ferri state and the PD state in an de- 
formed/ordered phase. If one of three sublattice order parameters vanishes and the other two 
remain finite, the PD state is considered to be realized. If three sublattice order parameters 
take finite values, the Ferri state is considered to be realized. 

We adopt the nonequilibrium relaxation method^^"^^) to detect phase boundaries. This 
method utilizes a relaxation function of a physical quantity in Monte Carlo simulations. First, 
an ordered state is prepared as an initial state of simulations. We perform a simulation and 
obtain a relaxation function of a physical quantity. Another simulation is performed changing 
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Fig. 3. A comparison between the Metropolis updating and the heat-bath updating. The NER func- 
tions of the structure factors behave in the same manner while the correlation time is about four 
times smaller in the heat-bath updating. An initial state is (lattice, spin)=(Ferri, PD). 



the random number sequence and another relaxation function is obtained. We take an average 
of relaxation functions over these independent Monte Carlo simulations. This sample average 
is free from a systematic error due to averaging correlated data. If the relaxation function of an 
order parameter exhibits an exponential decay, the system is judged to be in the paramagnetic 
phase. If it exhibits a converging behavior to a finite value, the system is in the ordered phase. 
The relaxation function exhibits an algebraic (power-law) decay at the critical point. The 
method has been verified to be particularly efficient in slow-dynamic systems with frustration 
and randomness. 

In the present simulation both spin and lattice variables at each lattice point are updated 
simultaneously using the heat-bath probability. We have checked that the Metropolis updating 
algorithm and the heat-bath updating algorithm yield the same relaxation behavior but with 
different correlation time. A simulation of the heat-bath is about four times faster as shown 
in Fig. 3. 

The most important point in the NER method is to exclude out the finite-size effect from 
the raw relaxation function. The method is based upon taking the infinite-size limit first. If a 
relaxation function is affected by a finite-size effect, the relaxation behavior misleads us into 
thinking that the temperature is in the paramagnetic phase even though it is in the ordered 
phase. A lattice size of the present system is L x (L + 1) x lOL with L ranging from L = 24 
to 120. We have checked the size dependence as shown in Fig. 4. The figure exhibits that the 
size effect appears quite early in the Monte Carlo step. A system with L = 120 is confirmed to 
be free from the finite-size effect until one hundred steps in this figure. For every parameter 
point of K, and temperature we have checked the size effect by comparing data of L = 120 
and L = 96 and confirmed a time range that the data are free from the size effect. Results of 
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Fig. 4. Size dependences of the structure factor of the spin. Data of sizes smaller than i = 96 notably 
deviate from data of i = 120. An initial state is (lattice, spin)=(Ferri, Ferri). 



L = 120 are presented in this paper. 

Initial conditions of the present simulations are mostly the ground-state configuration: the 
lattice structure is PG^cm and the magnetic structure is ferrimagnetic. This combination is 
referred to as (lattice, spin) = (Ferri, Ferri). We have also carried out simulations with initial 
configurations of (lattice, spin) = (Ferri, PD), (Ferri, Para), (PD, PD) and (Para, Para) in 
order to exclude out a possibility of initial state dependence in the final conclusion. 

Computations are carried out by a PC cluster with 14 nodes consisting of Pentium-4 and 
Athlon-XP CPU. Total computation time is eighty days. 

4. Results 

4-1 The 1/3 structure factor 

The first issue is to make it clear whether the structural phase transition and the magnetic 
phase transition occur at the same temperature or not. The NER of the 1/3 structure factor 
identifies the temperature that a lattice deformation or a magnetic order appears. 

Figure 5 shows NER plots of the structure factors at k = 1.0. We have started the 
simulation with an initial configuration of (lattice, spin)=(Ferri, Ferri). This combination 
is realized in the ground state. The NER of the structure factor of the lattice exhibits an 
exponential decay for T > 11.25. It exhibits a converging behavior for T < 10.0. The structural 
phase transition is considered to occur at Tc = 10.625±0.625. The spin is still paramagnetic at 
this temperature. The NER of the structure factor of the spin shows that the magnetic phase 
transition occurs between T = 7.5 and T = 9.0 {Tc = 8.25 it 0.75). Therefore, two transitions 
occur at different temperatures. Since the structural transition temperature is higher than 
that of the magnetic transition, the lattice system is considered to dominate the whole system 
at K = 1.0. 

We have carried out the same analysis changing a value of k as listed in eq. (7) and found 



8/15 



J. Phys. Soc. Jpn. 



Full Paper 



10^ 



10' 



1 1.00 



10' 







10^ 



t>0 



10" 



lattice 



10^ 



10 



t 

spin 



T= 5.Q0 
T= 6.25 
r= 7.50 

T= 9.00 
T= 10.00 



\r= 12.50 ■•r= 11.25 
10^ 



\ \ \ "\ T= 9.00 

\ \ /^idw 

7-12.5d \ \ \ K= 1.00 



T= 5.00 
r= 6.25 



r= 7.50 



10 
t 



10' 



Fig. 5. NER plots of the 1/3 structure factors of the lattice and the spin at «; = 1.0. The temperature 
is as denoted aside each plot. The structural transition and the magnetic transition occur at 
different temperatures. 



that the simultaneous transition may occur at k ~ 1.6 as shown in Fig. 6. Both structure 
factors decay exponentially for T > 11.25 and converge to finite values for T < 11.00. The 
transition temperature is estimated as Tc = 11.125 ± 0.125. For k = 2.0 and 2.5 it is found 
that the magnetic transition occurs at a higher temperature than a structural transition 
temperature. (Figures are not shown in this paper. Detailed data have been reported in the 
master thesis of one of the authors (T. S.).^^)) Therefore, there must be a point near k = 1.6 
and T = 11.125 that the simultaneous transition occurs. 

4-2 Sublattice polarization/magnetization 

We have observed the NER of sublattice polarization and sublattice magnetization in 
order to identify a lattice structure and a magnetic structure at low temperatures. It should 
be noticed that a sublattice may change its role of taking magnetization of t, i or depending 
on samples of independent Monte Carlo simulations. For example, when we start a simulation 
from the paramagnetic state at a temperature that the PD state is realized, one sublattice may 
take any three states of t, J, and in each Monte Carlo simulation. Therefore, a simple average 
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Fig. 6. NER plots of the 1/3 structure factors of the lattice and the spin at k = 1.6. The simultaneous 
transition occurs near T = 11.125. 



of sublattice magnetization over the independent Monte Carlo simulations yields an average 
of magnetization over sublattices: magnetization is obtained. In order to distinguish the above 
ambiguity, we collect data of sublattice magnetization depending on their value. Sublattice 
magnetization whose value is the smallest among three is stored in m^j^. For example, a 
spin-down sublattice (j) in the PD state contributes to this variable. That of the largest value 
is stored in m^gx (t in PD); ^'^d the rest is stored in m^j^j (0 in PD). 

Figure 7 shows the NER of squared sublattice polarization/magnetization at k = 1.0 and 
T = 6.25. The initial state is (lattice, spin)=(Ferri, Ferri). Three sublattice polarization con- 
verge to finite values. Here, m^;^ corresponds to j-state and its squared value is largest among 
three, m^j^ and correspond to t-state and converge to the same value. A realized state 

is considered to be the lattice Ferri (PGacm). The NER of squared sublattice magnetization 
suggests that the magnetic structure is the PD state, m^;^ and m^gx respectively correspond 
to the j-state and the j-state and converge to a finite value, m^;^ vanishes exponentially. 
Therefore, it is concluded that (lattice, spin)=(Ferri, PD) state is realized at k = 1.0 and 
T = 6.25. We have verified this result by starting simulations with (Ferri, PD) state, (PD, 
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Fig. 7. NER plots of squared sublattice polarization (lattice) and squared sublattice magnetization 
(spin) at K = 1.0 and T = 6.25. The lattice Ferri state and the magnetic PD state are realized. 



PD) state and (Ferri, Para) state. Every relaxation plot has drawn out the same conelusion.^^^ 
Figure 8 shows the NER of polarization/magnetization at k = 1.6 just above and below 
the simultaneous transition temperature. All the relaxation funetions deeay exponentially at 
T = 11.25. This evidenee suggests that both lattice and spin take symmetric/paramagnetic 
configurations at this temperature. Contrary, the relaxation functions converge to finite values 
at T = 11.00. The (lattice, spin)=(Ferri, Ferri) state is considered to be realized. This is 
the ground state configuration. A direct transition from the high temperature phase to the 
ground-state phase has occurred. Intermediate phases vanish at this point. 

5. Discussion 

Putting all the NER results together we obtain a T-k, phase diagram of the present model 
as shown in Fig. 9. When there is no spin degree of freedom at k = 0, successive structural 
phase transitions occur from the symmetric phase {P6c/mmc) to the lattice PD phase (P3cl) 
and then to the lattice Ferri phase (PGacm). Contrary, when there is no lattice degree of 
freedom at k = oo, we have also confirmed that successive magnetic transitions occur from 
the paramagnetic phase to the PD phase and then to the ferrimagnetic phase. 
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Fig. 8. The NER of sublattice polarization/magnetization at k = 1.6. The reahzed state is (lattice, 
spin)=(Para, Para) at T = 11.25, while it is (lattice, spin) = (Ferri, Ferri) at T = 11.00. 



In case where both lattice and spin exist, a structural transition and a magnetic tran- 
sition occur. If an energy scale of the lattice is larger than that of the spin {k < 1.3), the 
structural transition occurs at a higher temperature than magnetic transition temperatures. 
At this transition temperature a lattice changes its structure directly from the symmetric one 
to the ground-state one. The intermediate phase does not appear. The successive magnetic 
transitions occur at lower temperatures. This interesting phenomenon can be explained by 
the following scenario. A key concept is relaxation of frustration. 

A system with a larger energy scale becomes a master and the other one becomes a slave. 
When the structural transition occurs at a higher temperature, the lattice is a master and the 
spin is a slave. The spin as a slave serves a master relaxing frustration of the lattice by taking 
the PD configuration. As explained in Fig. 2, the spin favors the PD configuration when the 
lattice takes the Ferri structure, and vice versa. Since the temperature is still high, the PD 
spin configuration cannot be a long-range order. However, the short-range order is sufficient 
to relax frustration. The lattice can take the ground-state configuration without experiencing 
the intermediate PD phase. 

A phase boundary of the structural transition does not depend on k when k < 1.3, and 
smoothly interpolate to the phase boundary between Para and PD at k = 0. It is speculated 
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that the structural PD phase suddenly vanishes by an introduction of an infinitesimal spin 
degree of freedom. When the temperature is lowered, the successive magnetic transitions 
occur. There is no slave relaxing frustration of the spin variable. Therefore, the intermediate 
phase appears. 

Roles of the lattice and the spin are completely exchanged when an energy scale of the spin 
part exceeds that of the lattice part for n > 1.6. The spin becomes a master and the lattice 
becomes a slave. The magnetic transition occurs at a higher temperature and the ground- 
state configuration (ferrimagnetic state) is realized by an assistance of the lattice system. The 
successive structural transitions occur at lower temperatures. 

The simultaneous magneto-dielectric transition occurs when two energy scales coincide. 
Either system relaxes frustration of the other system cooperatively. Both realize the ground- 
state configuration at the same temperature. Intermediate phases of both systems vanish. The 
simultaneous transition temperature corresponds to the real temperature as 

{Tc/\Ji\) X Jf/K = 11.125 X 2.5/1.6 ~ 17K. (12) 

It is almost half the experimental result of 37K. This underestimate is due to adopting the 
large reductions of exchange integrals: A J^^ ^ 2) = 0-^'^(^i2)- '^^^ exchange integrals can be 
reduced to 60% of the original value. We made the reduction large in order to observe the spin- 
lattice cooperative phenomenon easily. However, it decreases an energy scale of the spin part, 
and consequently decreases the magnetic transition temperature. There is also an ambiguity 
in an experimental estimate of jf. Therefore, we consider that our estimate of the transition 
temperature is not inconsistent with the experimental result. 

This scenario possibly explains the simultaneous transition in RbCoBr3. It is also spec- 
ulated that there is no further magnetic transition at lower temperatures. The ground-state 
magnetic state (possibly the ferrimagnetic state) is considered to be realized below the simulta- 
neous transition temperature. Recently, it is reported that a single-crystal neutron diffraction 
experiment has revealed that all three sublattice magnetization are ordered at 37K?^^ This 
experimental evidence supports our speculation of an appearance of the ferrimagnetic state 
at the simultaneous transition point. 

In this paper we have introduced a model Hamiltonian which explains magneto-dielectric 
phase transitions in ABXs-type compounds. It is supposed that any deformation of a lattice 
decreases the magnetic exchange interaction. The model is made as simple as possible just 
to understand the phenomenon qualitatively. Choices of the system parameters are not real- 
istic. Therefore, there remains much to improve in order to compare theoretical results with 
experimental results quantitatively. 

An essential point of this phenomenon is relaxation of frustration by introducing two 
systems: lattice and spin. Frustration finds a way to be relaxed by utilizing the other system. 
If it is relaxed, the intermediate phase like the PD phase vanishes. It can be considered that 
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Fig. 9. A phase diagram of tlie present model. Symbols discriminate phases as denoted in brackets. 



the intermediate phase is a product of frustration that has nowhere to go. 
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